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Nanophotonic structures with irregular symmetry, such as quasiperiodic plasmonic crystals, have gained an 
increasing amount of attention, in particular as potential candidates to enhance the absorption of solar cells 
in an angular insensitive fashion. To examine the photonic bandstructure of such systems that determines 
their optical properties, it is necessary to measure and model normal and oblique light interaction with 
plasmonic crystals. We determine the different propagation vectors and consider the interaction of all 
possible waveguide modes and particle plasmons in a 2D metallic photonic quasicrystal, in conjunction with 
the dispersion relations of a slab waveguide. Using a Fano model, we calculate the optical properties for 
normal and inclined light incidence. Comparing measurements of a quasiperiodic lattice to the modelled 
spectra for angle of incidence variation in both azimuthal and polar direction of the sample gives excellent 
agreement and confirms the predictive power of our model. 



Coupled systems play a major role in nanophotonics and plasmonics. Plenty of applications such as optical 
nanosensors 1 " 3 are based on Fano resonances 4 which are due to interactions of different optical modes. 
The influence of planar periodic structures on particle plasmon resonances was studied by Lamprecht 
et al. 5 and Auguie et al. 6 Subwavelength holes in a metal film have been investigated by Garcia- Vidal and co- 
workers 7 . All of the mentioned publications have only examined periodic structures. However, also planar 
aperiodic and disordered structures have been studied 8 " 11 . Even 3D photonic quasicrystals were investigated 12,13 . 
The coupling between particle plasmons and waveguide modes has been studied in periodic 14 " 16 as well as in 
disordered 17 " 19 systems. Additionally, nonlinear experiments have been performed 20 " 22 . These systems are used 
for intriguing applications such as enhancing the efficiency of solar cells 23 ' 24 as well as subwavelength focusing of 
light 25 . The latter application is based on a quasiperiodic structure. One of the key properties of such quasicrystal- 
line systems is its more isotropic optical behaviour in comparison with periodic structures 26 . Therefore, the ability 
to model normally incident as well as angle -dependent optical properties is crucial. However, Fourier-periodic 
methods such as scattering-matrix calculations cannot predict the spectra of quasiperiodic and disordered 
structures. Therefore, the presented model of a previous paper for 2D plasmonic gold nanodot arrays on a slab 
waveguide for normal incidence 27 has to be modified and expanded for oblique light incidence. Our model utilises 
the Fourier transform of the spatial lattice positions in combination with the dispersion of the slab waveguide. 
Additionally, elliptically shaped nanodots are taken into account. Experiments with different incident polarisa- 
tions, different inclination angles as well as different azimuthal angles are compared to the theoretical curves, and 
excellent agreement is obtained. With the model presented in this paper, predictions of the optical properties of 
structures with a quasiperiodic gold disk arrangement are possible, which is confirmed by experiments. 
Therefore, the optical properties of future applications of these structures can be modelled and optimized by 
using this approach. Such applications include for example plasmon-assisted coatings for solar cells which would 
make their absorption more isotropic and polarisation insensitive. 

Results 

Modelling and simulations. In the model as well as the experiment, a sample is used which consists of a 180-nm- 
thick Hf0 2 slab waveguide on a quartz substrate with plasmonic gold dots on top, arranged in a quasicrystalline 
Penrose tiling (see Fig. 1 and Methods for details). Light with the vector k incident on the sample is defined by the 
angle of incidence 8 as well as the azimuth angle cp [see Fig. 2(a)]. In comparison with normal light incidence the 
incident k vector additionally possesses components parallel to the xy plane. Since the electric field vector E a is 
always perpendicular to k, E a can exhibit a component in z direction depending on (p and on the polarisation angle 
a between the component E a jXy and the x axis of the sample. Therefore, the theoretical model described in Ref. 27 
has to be expanded. First, we start by calculating the 2D Fourier transform of the structure. This is carried out by 
setting a Dirac delta function at each coordinate of the nanodisks 17,27 with x n andy n for the n th disk leading to a 
function 
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Figure 1 | SEM picture. SEM picture of the Penrose tiling of the plasmonic 
gold nanodots on top of an Hf0 2 slab waveguide. 
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The + sign in Eq. (5) belongs to the case when 0° < a — ft < 180°, 
whereas the - sign describes the case when 180° < a — /? < 360°. C a is 
a factor that depends on the incident wave. Note that only C TE is 
dependent on C a>z due to the fact that E a>z is always on a plane 
perpendicular to the direction of propagation. For dielectric 
nanodisks these are the only waves propagating in the waveguide 
and, therefore, the total intensity of the waveguide modes is given by 



rtot,TE z 



s p (k)c 2 TE dp, 



(9) 



and then taking the 2D Fourier transform of f(x, y) by using the 
projection slice theorem (see Ref. 28 and 29): 



(2) 



Due to the structure design, namely a grating on top of a waveguide 
layer, it is possible to excite the waveguide modes. When a transverse 
electric (TE) or a transverse magnetic (TM) polarised waveguide 
mode propagates in the waveguide slab, its energy corresponds to a 
propagation constant fi p defined by 30 



Pp = Ky + 



(3) 



with g describing the location of the Fourier components. This 
indicates that, in order to calculate the angular propagation 
constants, it is possible to shift the 2D Fourier transform by the k 
component parallel to the xy plane, k xy , which is displayed as 
vectorially shifted circles in Fig. 2(b). However, each Fourier 
component as well as k^y are dependent on the energy of the 
incident light. Therefore, for a specific polar angle 3, the 
magnitude of the vector k xy is different for each Fourier 
component. In order to find the correct angular propagation 
constants, the TE and TM waveguide dispersion relations given in 
Refs. 31 and 32 have to be used. In these dispersion relations, which 
are plotted in Fig. 2(c) as black solid (TE) and red dashed (TM) 
curves, the absolute value of f$ p given by 

P p = y/ khin(S) 2 +g* + 2k sin(% cos(^ - J~) (4) 

is taken. This allows for calculation of fi p in dependence on the 
corresponding k value of the incident beam for each Fourier 
component with distance g to the centre of the Fourier transform 
at an angle p FT between the direction of this Fourier component and 
the k x axis. The location of the angular propagation constants can 
then be easily calculated with Eq. (3). The directions of the waveguide 
modes are given as vectors starting from the centre of the Fourier 
transform and pointing to the different angular propagation 
constants. In each of these directions at angle p a TE as well as a 
TM mode can be guided, where the TE mode (TM mode) is defined 
with its electric (magnetic) field vector on a plane perpendicular to 
the direction of propagation. In order to find the correct amount of 
TE/TM waves propagating in direction the amplitudes of incident 
light with polarisation E a have to be split into the components [see 
Fig. 2(d)] 



Cte = ± 



Ci xy sm\a-f}) + Cl z , 



Ctm = C a , xy co$(a - /?) 
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where Sp(k^) describes the amplitude and phase information in 
direction p. The positive k^ axis of Sp(k^) is used for 0° < < 
180° and the negative for 180° < 0 < 360°. 

However, for metallic nanodisks we need to consider additional 
features, namely the particle plasmons. Light with an incident polar- 
isation E a can excite particle plasmons along the principal axes u, v, 
and z [see Fig. 2(e)]. This leads to the following components: 



Cpiu = C^xyCu cos(a-y), 
C p i v = C a>xy C v sin(a-y), 



(11) 

(12) 
(13) 



The values C w , C v , and C z in Eqs. (11) - (13) are dependent on the 
eccentricity and the size of the particle. The three excited particle 
plasmons can be considered as independently polarised light, which 
is incident on the sample. These electric field vectors (E p i Ui Ep/ V , and 
Epiz) can then excite waveguide modes in the same manner as 
explained above. Figure 2(f) shows how each of the electric field 
vectors of the particle plasmons E p i u and E p i v splits up into a TE 
polarised as well as a TM polarised wave in direction leading to 
the components 



Cteu = - Cpi u sin(^-y),Grev = C piv cos(/J-y), 
Ctez — Cpl z , 

Ctmu = C p i u cos(j3 - y) ,C TMv = C p iv sin(^ - y) . 



(14) 



(15) 



As above in Eqs. (5) and (6) for the electric field vector E a>z , the 
particle plasmon in z direction can only excite a TE polarised wave- 
guide mode. Note the - sign for the C TEu component which is due to 
the fact that E T£w points in the opposite direction than E TEv 
[Fig. 2(f)]. The complete components C TE , a dd an d Ctm, add are 
obtained by adding up the different components of Eqs. (5) and 
(14) for a TE wave as well as those of Eqs. (6) and (15) for a TM 
wave in direction: 



C 



TE,add : 



-- Cte + Cteu + Ctev + Ct 



CTE,add — Ctm + Ctmu + Ct 



(16) 
(17) 



The total intensity of the waveguide modes for metallic nanodisks 
can now be calculated by using 



with 



p360 
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(18) 
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Figure 2 | Simulation model, (a) Illustration of the incident k-vector on the sample, (b) Angular Fourier transform (filled circles) shifted by k xy with 
respect to the Fourier transform for normal incidence (open circles), (c) Dispersion relations for vacuum (green short-dashed) and quartz (blue dash- 
dotted) as well as for TE (black solid) and TM waves (red dashed) propagating in a 180-nm-thick Hf0 2 layer on quartz, (d) TE and TM waveguide modes 
in direction ft for incident polarisation E a>xy . (e) Particle plasmon excitations along the particle main axes, (f) Additionally excited TE and TM waveguide 
modes in direction /? due to the particle plasmons. 



(19) 



As already explained above, the positive axis of Sp(k^) is used for 0° 
< P < 180° and the negative for 180° < P < 360°. 



Ftot,TE as well as F totjTM consist of amplitudes at specific k values. 
The corresponding energies are obtained by using the above-men- 
tioned waveguide mode dispersion relations displayed in Fig. 2(c) as 
black solid (TE wave) and red dashed (TM wave) curves. As reference 
also the dispersion relations of the surrounding media are plotted 
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(vacuum: green short-dashed, quartz: blue dash-dotted). The ampli- 
tudes Aj of the main waveguide mode resonances in Eqs. (18) and 
(19) at the energies Ej are used in order to describe the transmission 
amplitude t by a Fano model 27 ' 33 ' 34 



t = A bg e u 



A p iu cos 2 (a-y)r p/M e z 'V 



A ph sin 2 (a -y)Y ph e i ^ 



E — E piv + iTpiy 



(20) 



where the sum over j models the waveguide modes with spectral 
widths Tj and phase fy. The first term is due to the directly trans- 
mitted wave with amplitude A bg = 1 and phase (j) bg = 0.085 n, 
whereas the second and third terms describe the particle plasmons 
in u and v directions with amplitudes A p iJA p i v , spectral widths T p iJ 

T p i v = 0.15 eV, and phases (j> p i u j 4> p i v = — . Note that the particle 

plasmon in z direction is neglected since the height of the nano- 
disk is much smaller than the diameters in u and v direction. 
Therefore, the particle plasmon energy in z direction is in a dif- 
ferent energy range as the measured spectra and the value C z of 
Eq. (13) is much smaller than C u and C v . The amplitudes A pXu and 
A p i v as well as the value A are fitting parameters so that the 
particle plasmon amplitudes and the waveguide mode amplitudes 
possess the correct ratio. The phase of a waveguide mode not 
coupled to a particle plasmon is (j)oo = 0.47 n, which is the case 
when the energy of the waveguide mode is far away from the 
particle plasmon energy. For low energies the phase of the wave- 
guide mode converges to 0oo, whereas for large energies it con- 
verges to (j)oo + 27i which corresponds to If waveguide mode 
and particle plasmon possess the same energy, <j>j has a phase-shift 
of 7i with respect to the uncoupled phase <^oo. Therefore, a phase 
behaviour of 



= 2 arctan ( _ I cos (a — y) 



+ 2 arctan 



(E E p ly 



(21) 



sin 2 (a -y) + 0oo+7T 



is assumed due to the two particle plasmons at energies E piu and 
Epiy. The extinction spectra can then be obtained by calculating 
Ext = — ln(T) with the transmittance T = \t\ 2 . 

Experiments. The sample described in the Methods section was 
fabricated and then measured with a white light transmission 
setup 15 . The measurement principle is also described in the 
Methods section. The angle of incidence S was varied between 0° 
and 6° in steps of 1° both in x {(p = 0°) and in y direction {(p = 90°) of 
the sample. All measurements, shown as black curves in Figs. 3 and 4, 
were performed for an incident polarisation with a = 0° and a = 
— 90°, respectively. The curves on the bottom of Fig. 3(a) and (b) as 
well as those of Fig. 4(a) and (b) were measured at normal light 
incidence. The broad resonances visible in all spectra without the 
dips are the excited particle plasmons. Due to the fact that the short 
main axis of the gold disks is rotated by 40° around the sample x axis, 
the particle plasmon resonances along both main axes are excited for 
incident polarisations with a = 0° and a = —90°. However, by 
comparing Fig. 3 to Fig. 4, one recognises that the particle 
plasmons are located at slightly different energies. This means that 
either the plasmon along the short main axis (Fig. 3) or the plasmon 
along the long main axis (Fig. 4) is more pronounced. The sharp 
resonances at approximately 1.5 eV, 1.6 eV, 1.75 eV, and 1.85 eV 
are due to TE and TM waveguide modes propagating in the 
waveguide layer with a propagation constant equal to the inner 
ring of the Fourier components [see inner dashed circle in 
Fig. 2(b)] as well as those with a propagation constant equal to the 
outer ring of the Fourier components [see central dashed circle in 



<p = 0°, a = 0° q> = 90°, a = 0° 



i — i — i — i — i — i — i — i — i — i — i — i — I I — i — i — i — i — i — i — i — i — i — i — i — r 
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Figure 3 | Measurements with polarisation a = 0°. Measured (black) and modelled (red) curves for a Penrose tiling with edge length d = 530 nm. The 
angle of incidence was varied from 0° to 6° (from bottom to top) (a) along the sample x direction (<p = 0°) and (b) along the sample y direction (tp = 90°). 
The incident polarisation was 0°. The curves are shifted upward for clarity. 
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Figure 4 | Measurements with polarisation a = —90°. Measured (black) and modelled (red) curves for a Penrose tiling with edge length d = 530 nm. 
The angle of incidence was varied from 0° to 6° (from bottom to top) (a) along the sample x direction {(p = 0°) and (b) along the sample y direction {(p = 
90°). The incident polarisation was —90°. The curves are shifted upward for clarity. 



Fig. 2(b)] . According to the energy position of the waveguide modes 
compared to the particle plasmons and thus to the degree of coupling 
between these resonances, the waveguide modes either appear as 
sharp peaks (decoupled, see resonance at 1.5 eV) or as sharp dips 
(coupled, see resonance at 1.85 eV). By tilting the sample, the 
waveguide modes split up into different components 16 ' 30 due to the 
different absolute values of the propagation constants for oblique 
light incidence. The behaviour for varying the inclination angle in 
x direction is different than that in y direction as expected. 
Additionally, we observe that this behaviour is dependent on the 
incident polarisation. This polarisation dependence has two 
reasons. First, the waveguide modes directly excited by the incident 
light show the same resonances with the same amplitudes for all 
spectra with a = cp = 0°, 36°, 72°, 108°.... The location as well as 
the amplitudes of these waveguide modes differ from the spectra with 
a = q> = 18°, 54°, 90°, 126°... due to the 10-fold symmetry of the 
Penrose tiling. Second, due to the fact that there are also waveguide 
modes present in the system that are excited by the particle plasmons, 
the amplitudes of the waveguide modes are also different for the 
spectra with for example a = q> = 0° and a = q> = 36°. However, 
the differences due to the second reason are only present, when the 
gold dots are elliptically shaped. The red curves in Figs. 3 and 4 are 
the spectra calculated with the simulation model presented in 
subsection Modelling and Simulations. A measured curve for one 
polarisation at normal light incidence was taken for fitting the 
theoretical spectrum to it. The fitting parameters obtained there 
were A piu = 0.0735, A plv = 0.1129, and A = 0.0016 and were used 
for calculating all other spectra. By comparing the measured and 
modelled spectra of Fig. 3, we observe an excellent agreement for 
both x and;/ variation of the inclination angle. The behaviour of the 
measured spectra of incident polarisation angle a = — 90° (see Fig. 4) 
with the same fitting parameters is also well reproduced by the 
calculated spectra which confirms our model. Measured and 
modelled curves are shown here for incidence angles up to 6°. 
However, our model can easily predict the optical properties at 
larger angles. 



Discussion 

We have introduced a theoretical model for simulating the optical 
properties of 2D metallic photonic quasicrystals at normal and 
oblique light incidence for all polarisations. This model includes 
the Fourier transform of the structure shifted by the k vector com- 
ponent parallel to the structure surface and all possible waveguide 
modes as well as particle plasmon resonances. The energy spectrum 
is obtained by utilising the waveguide dispersion relations and a 
phenomenological model. We have measured a quasicrystalline plas- 
monic crystal sample by varying the angle of incidence as well as the 
light polarisation and obtained excellent agreement between the 
measured and the modelled spectra. Our model could find applica- 
tions in plasmonic quasicrystalline super- oscillatory lenses 35 as well 
as plasmonically enhanced solar cells using waveguide-plasmon geo- 
metries 23 ' 36 . 

Methods 

Fabrication. The quasiperiodic lattice of the measured sample consists of 25-nm- 
high gold nanodisks which were placed on the vertices of a Penrose tiling with an edge 
length of 530 nm by using electron beam lithography. This arrangement is depicted 
in Fig. 1. Quartz (n = 1.46) with a 180-nm-thick Hf0 2 waveguide layer on top served 
as substrate. The refractive index of Hf0 2 is dispersive in the relevant spectral range 
and can be described by the following equation: 



n = 1.905- 



0.021 



0.0003 



(22) 



with X as the wavelength of the incident light in fim. The gold disks are elliptically 
shaped (97 nm and 118 nm long main axis diameters) and rotated by an angle of 40° 
between the short main axis of the particle and the sample x axis. 



Measurements. The light of a white light source was polarised and then focused on 
the sample by using a microscope objective (Zeiss, A-Plan, lOx, 0.25). The polariser 
was mounted on a rotation stage in order to change the incident polarisation. The 
light behind the sample was collimated and then focused on the slit of a spectrometer 
(Acton SpectraPro 500i with a CCD camera and a grating with 150 g/mm). A pinhole 
was used in the collimated beam after the sample in order to reduce the aperture angle 
of the beam below 0.2° 15 . The sample was placed on a rotation stage in order to 
perform measurements with different incidence angles on the sample. 
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